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'— ' Abstract 

I> In a set of experiments, Couder et. al. demonstrate that an oscillating fluid bed 

^ may propagate a bouncing droplet through the guidance of the surface waves. We 

pg present a dynamical systems model, in the form of an iterative map, for a droplet on 

CN an oscillating bath. We examine the droplet bifurcation from bouncing to walking, and 

C^ prescribe general requirements for the surface wave to support stable walking states. 

^^ We show that in addition to walking, there is a region of large forcing that may support 

the chaotic bouncing of the droplet. Using the map, we then investigate the droplet 
trajectories for two different wave responses in a square (billiard ball) domain. We 
show that for waves which are quickly damped in space, the long time trajectories in a 
square domain are either non-periodic dense curves, or approach a quasiperiodic orbit. 
C^ In contrast, for waves which extend over many wavelengths, at low forcing, trajectories 

tend to approach an array of circular attracting sets. As the forcing increases, the 
attracting sets break down and the droplet travels throughout space. 
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1 Introduction 

In a recent experiment ^, Couder et. al. demonstrate that under sufficient conditions, a 
fluid droplet will bounce infinitely often on an oscillating bed of the same fluid. Moreover, 
upon increasing the acceleration of the fluid bed, further experiments [H El IIZ] show a 
transition of the droplet from a stable bouncing state to a self propagating one. Here, the 
motion of the droplet is determined through a coupling with the underlying fluid bed. Upon 
each bounce with the fluid bed, the droplet generates a wave fleld in the bath. As the droplet 
continues to bounce on the bed, the previously generated waves act to guide the trajectory 
of the particle. Following [U [T7] we refer to the combination of the moving droplet dressed 
with a local wave as a walker. 

The resulting coupling between the droplet motion and underlying wave fleld has led to 
a variety of experiments IH El [71 El El [101 IH CSl [131 IIHI lEj demonstrating that the droplet 
trajectories may exhibit features reminiscent of a wave fleld (see also [2] for a review). For 
instance, in the flrst of a series of experiments [6l [17] , the authors show that two droplets may 
orbit or scatter without direct contact, but rather through the mediation of the underlying 
fluid bath. In addition, [17] goes further to investigate the phase diagram for bouncing and 
walking droplets with different sizes and accelerations of the fluid bed, as well as providing 
the flrst steps towards a phenomenological model for the droplet wave system. In their model, 
the authors average over one droplet interaction to obtain an ODE describing droplet motion. 

The experiments have not only been limited to the behavior of single droplets, but also 
include many droplet systems. For instance, [TT] demonstrate that droplet pairs may form 
localized, orbiting bound states, and even small crystalline structures. Further work [8j shows 
that the crystalline structures may include many different Archimedean lattices, while in [7j , 
Eddi et al. examine the instabilities in periodic hexagonal and square arrays. 

Other experiments have also analyzed the trajectories of droplets in bounded domains. 
For instance, |1] demonstrate that when a droplet is sent through a single slit scattering 
experiment, the droplet path may propagate with an apparently random scattering angle. 
Upon repeating the experiment for many droplets, the data shows that the droplets prop- 
agate with a probability distribution statistically analogous to a scattered wave amplitude. 
Meanwhile, in |9j Eddi et al. examine the motion of droplets in a conflned billiard setting. 
In analogy with quantum tunneling, they show that under sufficient conditions, instead of 
reflecting off the boundary walls, the droplets may occasionally cross a dead zone region 
which does not support stable walking trajectories. The paper also experimentally shows 
that one may obtain stable quasiperiodic trajectories, or the appearance of ergodic trajec- 
tories for droplets in a billiard domain. Lastly, [12J examine the trajectory of bouncing 
droplets in a rotating fluid. In a qualitative analogy with Landau orbits, they show that at 
sufficient forcing of the fluid bed, the droplets move in circular orbits with quantized radii. 
The papers [H Il2| also provide a simplifled phenomenological model for the droplet motion 
and numerically reproduce the qualitative behavior of the droplet trajectories. One draw- 
back of the model is that the fluid wave behavior is imposed as an ansatz without a model 



for the evolution of the fluid bath. Subsequent theoretical improvements were included in 
[To] where the authors provide a more detailed description of the waves generated near the 
Faraday threshold. 

Although in some of the previously mentioned experiments, the droplet trajectories ap- 
pear to approach stable structures, such as quasiperiodic orbits, in other situations, the 
trajectories appear ergodic or must be understood in a statistical sense. The goal of this 
paper is to provide a first step towards describing the long time behavior of droplet trajecto- 
ries, with an eye towards understanding their statistical behavior. We do this by adopting a 
dynamical systems approach, and model the motion of the droplets using a discrete iterative 
map. In particular we show that a map retaining one past bounce of the droplets history is 
simple enough to yield analytic results, yet complex enough to reproduce a variety of droplet 
behaviors: including the bifurcation from bouncing to walking, as well as quasiperiodic orbits 
and the appearance of dense trajectories in bounded domains. 

In the first section of our paper, we introduce an iterative map model for the droplet 
trajectories, and examine the general properties required for a bifurcation from bouncing 
to walking in free space. In the subsequent section, we introduce a gravity-capillary model, 
similar to the one introduced in [TU] , for the underlying fluid bath and analytically examine 
the bouncing to walking bifurcation. Using the wave model, we also show the possibility for 
chaotic droplet motion. Lastly, we examine droplet trajectories in bounded, square domains 
using two different wave models. In the first case, we use the gravity-capillary model and 
show the existence of either dense trajectories or quasiperiodic orbits. In the second case, we 
use a wave ansatz analogous to [U [12] and show trajectories may bounce between different 
regions of the domain. 



2 The Iterative map 

In this paper we are interested in understanding the two dimensional dynamics of the fluid 
droplets, as they propagate around the oscillating bath. As a result, we work in two spatial 
dimensions and record only the two dimensional position of the droplet, thereby ignoring 
the vertical (bouncing) motion of the particle. Let f2 C R^ be the two dimensional domain 
of the fluid bath, and denote the continuous time position and velocity of the fluid droplet 
by y(i) G Q and v(t) respectively. In addition, we let /i(x, t;y) denote the surface height 
of the fluid bath at x. Since the wave field depends on the past history of the particle, we 
include y(T) (for r < t) as a functional parameter in the wave field. 

To derive a set of dynamic equations, we adopt a simplified phenomenological model 
for the contact interaction of the fluid bed and droplet. In doing so we do not resolve the 
microscopic interaction between the fluid bed and droplet, but rather assume a spontaneous 
interaction. Specifically, we make the following simplifying assumptions pertaining to the 
interaction: 



Al. The nondissipative forcing from the fluid bath on the particle is proportional to the 
wave field slope at the time and location of contact. 

A2. At each bounce, the particle provides a point source forcing to the fluid wave field. 

We note that the general assumptions (Al)-(A2) are also made in the phenomenological 
model taken in [H |T2]. Following assumption (Al), the equations of motion for the fluid 
droplet take the form 

y = V (1) 

V = -[FVMx,t;y)|^^y + 7v-],5p(t), (2) 

where F is the amplitude of the forcing on the particle, and 6p(t) = 6p(t + T) is the periodic 
Dirac delta function. To capture dissipation in the droplet-bath interaction, we have added 
the additional term, 7v~5p(t), where v~ is understood to be the velocity prior to impact. 
Mathematically, we take w~{t)5p{t) = \(t — e)6p(t) with e — )■ 0, to correctly define the product 
of a distribution 6{t) with a discontinuous function v. One should note that in the general 
case, the shape of the wave field /i(x, t;y), depends on the previous history of the particles 
position y(t). 

Since the forcing is periodic, we may integrate the equations of motion ([I])-(|2| over one 
period to obtain a discrete map. To compactify the notation, we first choose the period 
of bouncing (T) as the natural time scale, and let y„ = y{n + e) and v„ = \{n + e) 
with (e — )■ 0), denote the droplets position immediately following the nth bounce. Upon 
integrating equations Q-Q, we obtain the iterative map 

Yn+l = Yn + ^n (3) 

v„+i = (l-7)v„-FV/i(x,n + l;yn,y„_i,...)| (4) 

' ^ — y n+l 

In the special case when 7 = 0, and h = /i(x) is independent of the droplets history, then ([s])- 
Q, reduce to the standard map [3]. The standard map is the discrete analogue of a particle 
moving in a strobed Hamiltonian system, and consequently the map preserves phase-space 
volume. We note that the presence of a path memory, with or without dissipation, breaks the 
discrete analogue of Liouville's theorem, thereby allowing stable attractors in phase space. 

Throughout this paper, we make a further simplifying assumption regarding the dissipa- 
tion in the droplet-wave field system. Specifically, we consider the limit when the dissipative 
time scale is much larger than the strobe period T, namely 

A3. The wave field depends only on the most recent droplet bounce. Explicitly, the as- 
sumption yields, /i(x, t;y„,y„_i, . . .) ^ /i(x,t;y„) for n < t < (n + 1). 

The assumption (A3) is therefore most applicable far below the Faraday threshold of the 
fluid bath where the long time path memory effects on the droplet are weak. 



With these assumptions, the iterative map becomes 

yn.+l = yn + ^n (5) 

v„+i = (l-7)v„ + Fg(yn+i,yn) (6) 

g(yn+i,yn) = -V/i(x,n + l;y„)| (7) 

Thus far, we have not exphcitly stated a model for the evolution of the fluid bath. 
Hence at this point one could take a variety of models for /;,(x, t; y) that accurately capture 
various properties of the underlying fluid field. In this paper, we first examine some general 
requirements for the wave field to support the stable walking of a droplet. We then study 
two explicit models for /i(x, t;y) in a square domain, with the goal of characterizing the 
resulting particle trajectories in two physical limits. 



3 General requirements for stable walking 

In this section we examine a set of general requirements for the map (|5])-([7]) to yield stable 
walking solutions in free space {Vt = IR?). 

Firstly, a droplet impact at position y„ will generate a radially symmetric wave field 
about Yn- Secondly, we note that the translational symmetry in free space implies that a 
wave field generated by an impact at y„ will only depend on the difference (x — y„). Hence, 
/i(x, t; y„) has the general form 

/i(x,t;y„) = /i(|x-y„|,t). (8) 

Letting r = |x — y„|, we can introduce the function g[r) = —drh{r, 1) where h{r, 1) is the 
wave field generated by an impact at time t = and evaluated at one strobe period later 
(ie. t = 1). Physically, g{r) describes the radial forcing on the droplet in the iterative map 
(|5|-([7]). Using the fact that y„+i — y^ = v„, we then obtain a one dimensional, iterative 
map for the droplet velocity in the radial direction 

Vn+l = il-j)Vn + Fg{Vn). (9) 

Without loss of generality, we align the system with the x-axis (ie. v„ = f„x) and further 
assume the parameters < 7 < 1 and F > 0. Since h{r,t) is radially symmetric, it follows 
that h{r,t) is an even function of r. Therefore, (provided h{r,t) is regular at r = 0), g{v) is 
an odd function of v and g{0) = 0. Hence, (tol) always admits f = as a fixed point. 

For stable walking solutions, we require the existence of a nonzero stable fixed point. The 
following criteria guarantee such a point. Let f * > and satisfy the following propositions 

PI. Existence of a nonzero fixed point 

g{v*) > 0. (10) 



P2. Stability of the fixed point 



v*g'(v*) 2 



:iii 



9[v 



7 



Here item (10) implies taking a forcing F = 'yv*/g{v*) > yields the fixed point velocity 
V*. Physically, condition (10) guarantees that the wave propels the droplet forward at each 



interaction. Meanwhile, condition (11) guarantees the stability of the linearized map at 
V*. Note that both conditions are invariant under rescaling g{v) — )■ Xg{v) (for A > 0). 



V 



Lastly, conditions (lO)-(ll) only guarantee a stable walking solution. An additional criteria 
for the bifurcation from stable bouncing to walking is an instability at v = 



P3. Instability of the v = fixed point 



v*g'iv*) 



g[v* 



> 1. 



(12) 



In the following section we examine a model for the wave field and show that dispersion and 
mild dissipation yield a wave field which supports stable walking solutions. We also remark 



that another standard wave field that fails condition (10) is the linear Green's function 



solution to the Helmholtz equation. Such a model is proposed in [^, however, they obtain 
walking solutions through the summation of many past bounces. 



4 A model for the wave field 



In the following section, we introduce gravity-capillary waves fT5^ as a model for the fluid 
wave evolution. We then examine the bouncing to walking bifurcation in free space, followed 
by the droplet motion in a two-dimensional square domain. In addition, we demonstrate 
that the wave field with weak fluid dissipation (/i < 1) supports stable walking solutions, yet 
does not support walking solutions for the strongly damped case (/i ^ 1). 

In the linear theory of gravity-capillary waves, one assumes an irrotational velocity field 
and models the fluid with a velocity potential (/)(x, z,t) with x G fi. Here, z is the vertical 
direction and z = is aligned with the unperturbed surface height h = 0. It then follows that 
/i(x, t;y„) and 0(x, z,t) satisfy the linearized gravity-capillary equations. Namely, (j){'x,z,t) 
is a harmonic function which vanishes at z — )■ — oo 

A0 = for 2 < 

(h = for z — )■ — oo. 



(13) 



In addition, 0(x, z, t) is coupled to /i(x, t; y„) via the kinetic and dynamic boundary condi- 
tions at z = 



ht 



bt + gh--/\2h + 2-<p,. 



(5(x - y„)5(t - t„) 

P 



(14) 
(15) 



Here A2 is the 2D Laplacian in x, while A is the 3D Laplacian in (x, z). Meanwhile, the since 
point source forcing at time tn = n enters as a delta function, via assumption (A2) and (A3), 
to model the instantaneous interaction with the bed. We remark that the dissipative term in 



equation (15) is only an effective dissipation as the assumption of an irrotational fluid field 
breaks down in a small viscous boundary layer near the wave field surface [TOl E 



m. 



To nondimensionalize the equations, we again use the period of bouncing (T) as the time 
scale, and take a length scalqjset by the pure gravity waves as L = gT'^. Letting, x — )■ Lx, 
t — )■ Tt, we also let — > (^0)(^), h —^ {Lh){^) where the dimensionless factor {^) is 
included to simplifjjj the equations to 



k + h- B'^A2h + 2fi(f), 



-5(x - y„)(5(t - tr, 



(16) 
(17) 



Here /i = v/{pg'^T'^) is the dimensionless viscosity, while B = gpL'^/a = g^pT^/a is the ratio 
of capillary waves to gravity waves. Typical experimental values [T7j, are a = .0209Nm~^, 
p = .96510^ kgm~^, T = 50~^s yielding i? ~ 7, /x ~ .03. Since we retain only the last position 



of the particle in the forcing, we supplement equations (16)-(17) with rest conditions at t 



(j){-K,Z,tn) 

/i(x,t„;y„) 




0. 



(18) 
(19) 



We therefore take equations (16)-(17) with initial conditions (18)-(19) to describe the fluid 



evolution between bounce n and ra + 1. 



Here equations (16)-(17), along with the boundary conditions, completely define an it- 
erative map for the droplets position and velocity. In addition, the map depends on four 
parameters, the dissipation of the wave (/x), the dissipation of the droplet bounce (7), the 
shape of the dispersion relation {B) and the acceleration or forcing of the particle {F). In 



the subsequent sections, we study (16)-(17) to understand the resulting particle motion for 
various forcing and domains. 



4.1 Solution in free space 

In this section we construct the single bounce iterative map (pi) in free space {Vt = R^). We 
do this by solving for the wave field /i(x, t;y„). With the exception of adding dissipation, 
the solution follows very closely to that of the standard Cauchy-Poisson problem for forced 

*We remark that another natural (Faraday) length scale is set using the linear dispersion relation w^ = 
gk + ak^ / p in combination with the period T. Such a choice, however leads to a more complicated form for 
the coefficients in the dimensionless equation without any difference in the underlying physics. 

^We may also reabsorb the factor ('^^) into the forcing amplitude in mh by nondimensionalizing F — > 



gravity-capillary waves. Specifically, we seek an expansion for of the form {k = |k|) 

(/.(x, z,t)= f A(k, t)e^'e'^'' dk. (20) 



Upon substituting the ansatz (20) into equations (16)-(17), one obtains an initial value 



problem for each y4(k, t). The solution for h{'x,t;yn) over t„ < t < tn+i then becomes 

/z(x,t;y„) = - /-^sin(a;B(t-t„))e^'^("-''"^e-^''(*-*"Mk, (21) 

J Ztcud 

where cojj = {k + B~^k^) — ij?k^ is the dispersion relatioirl To obtain an expression for the 
iterative map, we evaluate /i(x, t; y„) at time t^+i = t„ + 1 to obtain 

/i(x,t„+i;y„) = - /'-^sin(a;B)e*''(--y")e-^'='dk. (22) 

J 2tcud 

We now remark that one may asymptotically approximate /i(x, t„.+i; y„,) for yU ^ 1. The 
result is an overdamped wave, which does not support the steady walking of droplets. For 
large /i, the terms inside the integral may be approximated as follows: 

/i(x,t„+i;y„) ^ _y_l_e-iV'^(-y")dk + 0(/i-3/2) (23) 

47r/iyo Jo 

' +0(^-3/2). (25) 



27r Jq {B — t2jj,r cos i 

Here, we have aligned the droplet position and velocity with the x-axis (y„ = r„x, v„ = t;„x) 
and introduced r = |x — y„|. To compute the last integral, we make the change of variables 
z = e^^, and proceed by evaluating the residues enclosed by the unit circle \z\ = 1. For large 
fj, we have the over damped wave: 

%,Wi) = - ;^ ,^,, +0(/.-^/^). (26) 

y/B-^ + 4(/ir)2 



We also note that the error term 0{fi '^''^) < 0.5/i ^'^, so that (26) provides a good approx- 
imation for values of fi as low as /i ~ 5. Taking drh{r), we have 

n ^ 4/i2Ft;n 

"'^^^ = ('-^)""- (5- + 4^.2.^)3/2 • ^'^^ 



The iterative map (27) has only one fixed point Vn = 0, regardless of F and B. Hence, the 



strongly overdamped waves do not produce a bifurcation to walking motion as one increases 



411 the case when ujy < 0, the function 5iBi£sig ^fe becomes a strict exponential decay, corresponding 
to overdamping of the large modes. 



F. Physically, the effect of damping smooths out the wave curvature, thereby inhibiting a 
transition to stable walking. 

We now examine the wave field /i(r, t„+i) for small values of yU. To compute the field, we 



first integrate over the angular component 6' in (22) to obtain 

r-oo ^2 



/i(r, t„. 



+1) 



^D 



sin{ujD)Jo{kr)e ^'' dk 



{2t 



As a one-dimensional integral, we may numerically evaluate (28) and compute the iterative 
map 



Wn+i = (1 - j)Vn + Fg{Vn) 

g{v) = / — sm{ujD)Jo{kv)e-f'^ dk 

Jo ^D 



(29) 
(30) 



where Jq{z) = -^Joiz)- We note that f„ = is a fixed point of the map (29) for every value 
of B and fi. Hence, bouncing droplets with a fixed location r^+i 

of pgl). 



Tn are always solutions 



To illustrate the nature of the bifurcation from stable bouncing to walking exhibited in 
(29), we fix values of 7, /i and B, and continually vary F as the bifurcation parameter. As 
an explicit example, we take B = 7 in analogy with experiment, and /x = 1/4000 to be 
very smalPl With these parameters, one may then approximate the wave field h{r, t) with a 



polynomial to model the dynamics. 

Vn+l -- 

9{v) -- 



(1 - -f)Vn + Fg{Vn) 

llw(l + .76t;^-2.18t;^) 



(31) 
(32) 



< 



The polynomial fit (32) accurately reproduces the wave field over velocity values |f„ 
1.2. Here figures (la) and (lb) show the wave field and Poincare map, with the associated 



polynomial fit ( 32 ) . 



As one increases F, the fixed point solution f „ = becomes unstable at which point the 
system undergoes a pitchfork bifurcation. For values of F above the critical forcing {F > Fc), 
the points Vn converge to a new fixed point solution v*, indicating a bifurcation from stable 
bouncing to walking. For instance figure (2a) shows the bifurcation diagram in the case 
when 7 = 1. As one further increases F, the stable walking solution bifurcates a second time 



into a two-period orbit, followed by a transition to chaos (figure (2b)). We remark that the 
pitchfork bifurcation and transition to chaos occur for a large range of fi and B, however the 
subcritical bifurcation seen in (2a) is not generic. For instance, by varying B (ie. the shape 
of the wave) one may realize both positive and negative values of drrrrh{0), the relative sign 
of which determines whether the bifurcation is sub or supercritical. 



^The wave field h{r, t) for 7' < 1 is only weakly depended on yu for small values of /i. 





(a) Gravity-capillary wave field. (b) Poincare map. 

Figure 1: Shows the wave field and Poincare map for gravity-capillary model {B = 7, fi '^ 1] 



The dashed line shows the polynomial fit (32). 
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(a) Bifurcation diagram. 
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(b) Droplet velocities versus forcing. 



Figure 2: Shows the bifurcation from stable bouncing to walking for gravity-capillary model. 
In (a), the transition is a subcritical pitchfork where Fc = 11, yU ^ 1. In (b), after the initial 
bifurcation, the droplet undergoes a transition to chaos. 



4.2 Movement in a square 

In this section we examine the movement of droplets in a square domain using the model ([s])- 
(|6| with the wave field described by equations ( 16 )-( 17). To capture the refiection of the fiuid 
waves against the wall, we impose a Neumann boundary condition on the velocity potential 
^ = corresponding to a no fiuid fiux boundary condition on the bath. Equivalently, such 



a condition corresponds to a Neumann boundary condition on /i(x, t; y„): 

dh ^^ , ^ 

— = 0, xedQ 33 

an 

where n is the unit normal along the domain boundary. For instance, differentiating equation 
(16) and projecting onto the boundary yields 9t^ = dz-^ = 0. Hence the boundary condition 



j^= implies ^ = C, a constant in time. Since the constant C = at t = 0, we take 
1^ = for all time. 

an 

To evaluate the motion of the droplet in a square domain, we must compute the wave 
field /i(x; y„) at each iteration of the map. To aid in the computation of the field, we may 
exploit the method of images p] and the geometry of a square. For instance, since the field 
/i(x; y„) satisfies Neumann boundary conditions, and is generated by an impulse at each 
step, the solution may be generated by an infinite array of image points of the free space 
wave field hpsi^) at properly chosen locations x*™. The location of the image points depend 
on y„, the closest ones being at points reflected across the domain wall boundaries. One may 
then compute the wave field from the knowledge of the radial free space impulse response 

roc 7,2 

hpsir) = - / —sm{uD)Mkr)e-'''' dk (34) 

h{^;y^) = /i^<,(|x-y„|) + ^/i^5(|x-x}-|). (35) 

j 

Figure ^ illustrates the wave field with the most important image points. The addition of 
the image points yield correct boundary conditions for /i(x; y„) on the bottom and right side 



of the square. Since hpsi^) typically decays quickly, one may truncate the sum (35) for the 
efficient computation of the wave field. In our case we keep the first order contributions as 
illustrated in figure ^. 

In our numerical evaluation of the map, we fix 7 = 1, the dispersion B = 7 and the 
(dimensionless) size of the box D = 10 to be much larger than one wavelength. We then 
examine trajectories for different fluid dissipation < /x < 1 and forcing F. Since we model 
the droplet interaction with the boundary of the domain entirely by reflected waves, at large 
forcing for different parameters /i, there is a possibility that the droplet may physically 
collide or jump over the boundary. We therefore limit our attention to parameters which 
yield bounded trajectories, namely those which reflect off the walls. For example, figure ^ 
shows the allowable region of parameters to yield bounded walking solutions for different 
values of /x. 

Depending on the parameters of the underlying wave field, the long time behavior of 
particle trajectories may be classified into three categories depending on the nature of the 
limiting set: those which approach a flat set, those which approach a circular quasiperiodic 
orbit, and those which continually traverse the domain. In the last category, the trajectories 
appear to form a dense set throughout the spatial domain. To illustrate the different sce- 



narios, figure (5a) shows a trajectory collapsing towards a flat set where the droplet bounces 

10 



between two opposite walls, while figure (5b) shows the path of a droplet approaching a 
circular quasiperiodic orbit. As an example of the last case, figure (5c) shows part of a 



dense trajectory, while (5d) is the corresponding time averaged density plot. The dark ar- 



eas indicate positions the droplet visits more frequently than light regions. The emergent 
pattern, however, is not related to cavity modes of the square, but rather results because 
droplet trajectories tend to travel along paths near angles of 7r/4 with respect to the x-axis. 
For instance, although not shown, the statistics of the droplet velocity angles are centered 
around angles of ±7r/4. 




Figure 3: Plot shows the closest image points used to compute the wave field in a square 
domain. Each image point acts as a source with wave field hpsi'f')- 



5 A Long range model wave field 

As one increases the forcing of the fiuid bed towards the Faraday threshold [S], impulses 
made on the fluid surface generate waves with a fixed wavelength. In this section, we 
examine the long time behavior of droplet trajectories for an idealized response of the fiuid 
bed with a fixed wavelength. We note that we still adopt assumption (A3) to provide a simple 
analytical model, however, as one approaches the Faraday threshold, the wave field depends 
on a significant number of past interactions. These interactions can lead to additional effects 
on the droplet trajectory. 

In analogy with [H [12] we provide a phenomenological ansatz for the radially symmetric. 



11 




0.25 



Figure 4: Parameters for the gravity-capillary wave model in a square domain [D = 10). 
The solid line shows the threshold for walking droplets, while the dashed line shows the 
critical forcing for trajectories to remain bounded. The squares show two parameter values 



which yield flat attracting sets (fig. (5a)), the circle corresponds to a circular quasiperiodic 



orbit (fig. (5b)), while the triangle corresponds to a dense trajectory (fig. (5d)). 



free space, wave impulse response which decays with a fixed wavelength. Here we take 



hpsir) = Jo{kfr)e 



-fir 



(36) 



which differs from [U [12] by remaining regular at the origin (r = 0). As a result, for regions 



of the free parameters, the wave (36) satisfies the general requirements for droplet walking. 



To compute the wave field in a square domain, we again use (35) with the field (36), and 
propagate the droplet trajectories according to the map ([5])-([7|) taking 7 = 1. Here kf is 
the Faraday wavelength and fi is the spatial dissipation. In this section we take kj^ as the 
natural length scale so that the variables rescale as r — )■ kj^r, /i — )■ kj/j and D — )■ kJ^D is 
the dimensionless box size. In the following we also that 7 = 1, so that the iterative map 
has three free parameters: dissipation /i, D the domain size relative to the Faraday length, 
and F the forcing on the droplet. The dimensionless wave field is then 



hpsir) = Mr)e 



-fir 



(37) 



To determine the long time behavior of the map, we fix a set of parameter values and 
examine the trajectories for many different initial conditions. The data for each initial 
condition is chosen to survey the phase space within a bounded set by prescribing, |vo| < 0.25, 
and taking a maximum distance between yg and the nearest wall to be less than 1. Here the 
exact bounds of 0.25 and 1 are chosen somewhat arbitrarily to include a large, physically 
relevant, region of phase space, while avoiding unphysical difficulties (such as unconfined 
solutions). 
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(a) 7 = 1, F = .36 and ^ = 1/12 
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Figure 5: The long time spatial trajectories for section 4.2 collapse into flat (a), and circular 
(b), quasiperiodic orbits. Meanwhile trajectories in (c) travel throughout the domain. Figure 
(d) shows the long time spatial density (total time t = 10^) for a trajectory from (c). 



13 





(a) F = 0.5 



(b) F = 1.0 




(c)F 




1.1 



(d) F = 1.2 





(f) F = 1.4 



Figure 6: Attracting sets for the model response (36) and various forcing. Here D 
/x = l/10. 



8 and 



For each initial condition, we remove any transient effects by first evolving the trajectory 
for several thousand iterations. After discarding the transients, we then evolve the droplet for 
several thousand more iterations, and project the trajectory from the four dimensional phase 
space (x, v) onto the two dimensional physical domain (x). In many cases the trajectory 
approaches an attracting set in the form of a periodic or quasiperiodic orbit. We now describe 
in detail the long time trajectories as one varies F for fixed /i and D. Specifically, we consider 
in detail the case of Z) = 8, which corresponds to a box size of roughly one wavelength, 
and D = 12 which is just under two wavelengths of the wave field. In both cases, we fix 
/i = 1/10. Although varying the dissipation may alter exactly when bifurcations occur, or 
even the number of stable attractors in the dynamical system, it does not tend to change 
the qualitative behavior of the attracting sets. 

Initially, with small values of the forcing F, the long time trajectories approach one of 
several quasiperiodic orbits. These orbits form attracting sets for different regions of phase 



14 



space. For instance, the exact orbit a trajectory approaches depends only on the trajectories 
initial conditions. Together the collection of all attracting orbits form a symmetric array on 
which the precise pattern depends on the parameters F, fi and the box size D. Qualitatively 
the number of orbits depends most strongly on the box size D. The reason is that the droplets 
tend to localize near the troughs from the waves reflected off the domain boundaries. The 
larger box sizes allow for more wavelengths from the reflected waves. For instance, with 
fi = 1/10 and a wide range of forcing F, a box size oi D = 8 supports 4 quasiperiodic 
attracting regions, while D = 12 contains 8. 

As the forcing increasing, the spatial radii of the quasiperiodic orbits grow. For instance 
figures ^ and (JTl) show the attracting sets for box sizes 8 and 12 with different forcing. At 
sufficient forcing, the nature of the attracting sets change from thin circular orbits to thick, 
sets. At large forcing, the localized attracting sets break down, and the droplet wonders 
throughout the domain. Figures ([6JF) and ([7]f) show a shaded probability distribution for 
the droplets position at the large forcing. Despite the fact that there are no longer circular 
quasiperiodic orbits, the droplet still spends a significant time near the former quasiperiodic 
orbit regions. For instance, there are similarities in the dark outlines of ([6fe) and ([6]F), as 
well as (iTJl) and (ItIF). Although the distribution shows dark, highly traversed regions, and 
light, vacated regions in a regular array that appears related to the underlying wave field, the 
exact dependence is not completely understood. Here we defer further investigation to future 
work. For instance, in future work we seek to examine the relation between the dynamical 
systems invariant measure, and the underlying wave field. 



5.1 Conclusions 

Through the introduction of an iterative map, we model the dynamics and trajectories of 
bouncing droplets on an oscillating fiuid bed. As a first step, we examine the droplets 
bifurcation from a stable bouncing state to a stable walking one. In addition, we list several 
requirements for the underlying wave field to support stable walking droplets. Using the 
map, we then investigate the droplet trajectories for two different wave responses in a square 
(billiard ball) domain. The two responses are set to capture two different physical regimes: 
one where the wave field decays quickly, and the other where it extends over longer spatial 
length scales. In the first model, we provide a detailed analysis of the bouncing to walking 
bifurcation. We also show that, provided the droplet does not coalesce with the bed, there is 
a region of large forcing that may support the chaotic bouncing of the droplet. In the case of 
the square domain, we recover limit cycle and space filling trajectories which appear similar 
to those reported in |3]. Lastly, in the second model we show that for low forcing, trajectories 
tend to approach circular attracting sets. As one increases the forcing, the attracting sets 
break down and the droplet tends to travel through space, jumping between the former 
attracting regions. In future work we plan to further examine the statistical nature of the 
droplet trajectories, including their transport properties and invariant measures. 
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Figure 7: Attracting sets for the model response (36) and various forcing. Here D = 12 and 
/x = l/10. 
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